hep-th/9510131 
SINP-TNP/95-15 



Exact Bethe ansatz solution of nonultralocal quantum 

mKdV model 



Anjan Kundu 

Theory Group, Saha Institute of Nuclear Physics, 
1/AF Bidhan Nagar, Calcutta 700 064, India, 
e-mail: anjan@saha.ernet.in 



Abstract 

A lattice regularized Lax operator for the nonultralocal modified Korteweg de Vries 
(niKdV) equation is proposed at the quantum level with the basic operators satisfying a 
g-deformed braided algebra. Finding further the associated quantum R and Z-matrices 
the exact integrability of the model is proved through the braided quantum Yang-Baxter 
equation, a suitably generalized equation for the nonultralocal models. Using the alge- 
braic Bethe ansatz the eigenvalue problem of the quantum mKdV model is exactly solved 
and its connection with the spin-^ XXZ chain is established, facilitating the investigation 
of the corresponding conformal properties. 



1 Introduction 



The theory of quantum integrable systems made a major breakthrough with the formulation 
of quantum inverse scattering method (QISM) and the algebraic Bethe ansatz These 
techniques enabled us to solve exactly the eigenvalue problem of a number of quantum models 
including the spin chains p[ as well as the field theoretic models 0- 0. However, this success 
apart from some recent developments was confined mostly to a limited class of models, known 
as ultralocal models. For such systems the Lax operators at different lattice points j, k]jy^k 
must commute, i.e. 

Lbk{lJ')Laj{X) = Laj{X)Lbk{fJ-), aj^b, (1.1) 

where a, b signify associative space indices. Due to this constraint the basic equation given by 
the quantum Yang-Baxter equation (QYBE) 

Rab{X - fl)Laj{X)Lbj{fl) = Lbj{fi)Laj{X)Rab{X - /i) , j = l,...,N (1.2) 

can be multiplied for different j's and raised to its global form 

Rab{X - fl)Ta{X)Tb{fl) = Tb{tl)Ta{X)Rab{X - /i) (1.3) 

through the monodromy matrix 

Ta{X) = LaN{X)La,N-l{X) . . . L„i(A). (1.4) 



Taking trace of ( |1.3|) one derives in turn the relation [trT(A), trT(/i)] = 0, and since ln(trT(A)) = 
Yl'^=i generates the conserved quantities, one gets [C„, Cm] = 0, for n ^ m. That is the 
infinite number of conserved quantities {C^} are in involuion, which ensures the integrability 
of the corresponding quantum system represented by the Lax operator Lj. 

However, since the above formulation of quantum integrability depends crucially on the 
ultralocality condition ( |1.1|) , it fails for the nonultralocal models like WZWN, nonlinear a- 
model, Korteweg de Vies (KdV), modified KdV (mKdV) models etc. not satisfying (|1.1|) . 
The Lax operators Lj of such models involve discrete fields not commuting at different lattice 
points, which at the continuuum limit may correspond to commutation relaions like 

[v{x),v{y)] = 6'{x - y) 

expressed through derivatives of the 5-function 0. 

Nevertheless there were some progress in treating a few of such nonultralocal models [Q- 
131 as well as in formulating the basic theory |]T^- |T^, though the situation in dealing with 



genuine quantum integrable models is still far from its total perfection. In this context solving 
quantum KdV and mKdV models at par with the elegant approach of QISM |jl|, ^ remained 
as a major challenge. 

Our aim here precisely is to address this problem by presenting a quantum mKdV model 
and solving its eigenvalue equation through algebraic Bethe ansatz, much in parallel to the well 
established treatment of sine-Gordon ||^ or Liuoville 0] model. It should be mentioned that 
in a very recent work an indigenous nonultralocal formulation of KdV equation as a quantum 



mapping was given by Volkov [18 
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2 Classical mKdV model 



Before presenting the quantum mKdV model let us look briefly into some of its classical aspects 
as a well known integrable system [19| given by the field equation 

± Vt{t, X) + Vxxx{t, x) - 6v'^{t, x)vx{t, x) = 0, (2.1) 

and possessing an infinite set of conserved quantities C„, = 1, 2, . . . including the Hamiltonian 

H = Cs = -J_Jx{vl + v'). (2.2) 

The Poisson bracket (PB) structure for the field is given in a noncanonical form 

{v{x),v{y)} = T6'{x-y), (2.3) 
using which equation (|2.1|) can be derived from (|2.2|) as a Hamilton equation. Different signs 



in equation (2J) corresponds to fields with different signs in their PB relations (2^). Note 
that this noncanonical Poisson bracket (|2.3|) is the source of the nonultralocal behavior of the 
mKdV model at the quantum level. 

As in all classical integrable systems the nonlinear equation ( p.l|) can also be obtained from 
the linear system $3. = [/$, $f = as a flatness condition 



with the Lax pair 



Ut-Vx+[U,V]=0 



U{v,0 = ^(_l (2.4) 



and 



/ Ate + 2Cv^ - + 2iivx + ivxx - 2z^;3 \ 



It should be remembered that among this pair the Lax operator U{v,^) plays more important 
role and can generate alone all conserved quantities, which with the subsequent application 
of the PB derives the given equation from the Hamiltonian and other hierachies from the 
higher conserved quantities. The PB structure associated with the Lax operator also carries 
important information about the classical r-matrix. 

There exists an important connection between mKdV (|2.1| ) and the KdV equation 



± ut{t, x) + Uxxxit, x) - 6u{t, x)ux{t, x) = 0, (2.6) 

through Miura transformation u = Vx + v"^ which maps the PB relation (|2.3| ) for u{x) = 
2^ J2n ^nC"™^ to the classical Virasoro algebra 

H^n; ^m} = ('^ ~ ^Yn+m + ^f^n+m.O (2-7) 
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3 Lattice regularized quantum Lax operator 



We will be interested in the quantum generalization of the mKdV model. However our point 
of focus will be not the operator evolution equation, but the algebraic aspect of the problem 
related to its eaxct quantum integrability. For this purpose we consider a lattice regularized 
Lax operator 

involving basic quantum operators W^, which satisfy an interesting g-deformed algebra 

W+W- = qW-W+, [^t^t^ = 0' (3-2) 



with q = e*^''' , similar to the deformed braided algebra . For checking that at the continuum 



limit one recovers the relations relevant to the mKdV fields, we introduce operators Wj = e™j 
with the commutators ^ 

bfc > v^] = Til-{Sk-i,i - Sk,i-i), (3.3) 



and 



h 

] = - 25fc,/ + 6k,i-i), (3.4) 



which are consistent with algebra (|3.2| ). 

At lattice constant A ^ 0, the discrete operators go to the quantum field Av^{x) 
and due to the limit ^xix — y) the relation ( |3.3|) reduces clearly to the nonultralocal 

commutation relation 

[v^{x),v^{y)] = ±ij^{S^{x -y)- 5y{x - y)) = ±ijhS'{x - y) (3.5) 

recovering at the classical limit the PB relation ( p.3|) with ± signs corresponding to the fields 
v^{x). The reduction of operators ~ l + iAv^{x) at the continuum limit yields from ( |3.1| ) 

with 

C{v-,i) = -i{v-{x)a^ + ia^) (3.6) 

linked to the well known mkdv Lax operator (|2.4] ) as U (f , ^) = a^C{v~ , ^) involving field v{x), 
which is related to equation ( p.l|) with + and PB (|3.3|) with — signs. It is remarkable that 
though the lattice Lax operator ( |3.1| ) contains both variables having nontrivial commutation 
relation ( |3.4| ), the dependence on one of them drops out from the continuum Lax operator (in 



the present case it is v~^). Moreover since ( p. 4] ) at A — yields [v^ (x) , v~ (y)] = 
A ■ dxx^ix — y) — > 0, the fields answering to different signs of mKdV equation and the PB 
relations get decouppled at the continuum limit. Note that interchanging the operators 



3 



W along with q ^ q ^, keeps the underlying algebra (|3.2| ) invariant, while the corresponding 
Lax operator is changed as Lk{v~ ,v~^ , ^) Lk{v~^ ,v~ ,^), with the continuum limit 



where C{v^,^) dependents now on the complementary field v'^ related to the — sign of (|2.1|) . 
We suspect that this feature is likely to be present in the lattice regularization of all such 
nonultralocal nonrelativistic models exhibiting g-deformed algebra. 

4 Braided QYBE and quantum integrability of mKdV 
model 



For establishing the quantum integrability of the model associated with the Lax operator ( |3.1| ) 



we have to generalize the QYBE (|1.2|) to its braided variant ||T5l 

RatiX - fi)Z^^'L,j{X)L,,{f,) = Z-,'L,,{f,)La,{X)Rab{X - /^), J = 1,...,N (4.1) 
complemented at the same time by the braiding relation 

Lbj{lj)Laj+liX) = Laj+l{X)Z-,^Ujifi), [Ujifi),La^{X)] =0, t > J + 1, (4.2) 

which refiects the nonultralocal nature due to the noncommutation of the operators at neigh- 
boring lattice points. In such nonultralocal models along with the standard i?-matrix solution 
of the Yang-Baxter equation (YBE) 

Rl2{X - //)/?i3(A)i?23(/i) = i?23(/^)i?13(A)i?l2(A - fi) (4.3) 

one needs to find a Z-matrix satisfying similar equations 

i?2l(A)-Z'i3Z23 = Z232'i3i?2l(A) Zi2Zi'iR^2{X) = i?32 (A) Z'i3Z'i2 (4.4) 

and 

■Z'l2^13^23 = Z23Z13Z12. (4-5) 

In general the Z-matrix may also depend on the spectral parameters. Note that the set of 
YBE's ( t4.3D -( ^75D is compatible with the associativity condition of ([4.1|) . This formulation is 
specialized to handle the nonultralocality not extending beyond the next neighbouring points, 
which agrees with the present mKdV model with the field commutator ( p. 51 ). 
We may define the monodromy matrix in the interval as usual by 

ri'='^](A) =L„fc(A)L„,fc_i(A)...L,,(A) 

acting in the space 7i = ® Vk-i ® . . . V,-. It is important to note that due to the braided 
QYBE ([4.1|) and braiding relations ([4.2|) the transition to the global equation (in analogy with 



the transition from (]_^) to ( p..3|) ) is possible also for this nonultralocal case resulting 
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(4.6) 

where the similarity in form with (|4.1| ) reflects the underlying Hopf algebra structure. 

For physical models with periodic boundary condition Laj+N^X) = Laj{X) more interesting 
thing happens to the global monodromy matrix Ta{X) = T]^'^1(A) for the closed chain [A^, 1], 
since now Li{X) and Ljv(A) becomes nearest neighbours again with nontrivial commutation 
relation (^.2|). The corresponding global braided QYBE is given consequently by 



Ru{X - fi)Z^,'T,{X)Z{,'UiJ,) = Z^^^TMZ2iT,{X)Rr2{X - fi) 



(4.7) 



Equipped with this theory for quantum nonultralocal models we may start now with the lattice 
Lax operator ( p.l| ) of the mKdV model, which from the braiding relation ( |4.2| ) and with the 
use of commutators ( |3.2| ) for neighboring points yields a solution 



n2 



^21 



For finding the quantum i?-matrix we plug Z-matrix 
QYBE (gl|) and use relation like 



(4. 



together with in the braided 



{W^)''iWff = g"^(M^-)'^(iy+)" 
for the same lattice point along with the trivial one 

After some algebra we arrive at the remarkable result that the i?-matrix of the lattice 
mKdV model is given by the same trigonometric solution 



R{X-fi) 



( a(A — /i) \ 

6(A — y) c(A — jj) 

c(A — ii) b{X — fi) 
V a{X-y) J 



(4.9) 



with 



a(A — fi) = sin(A — /i + ^7), b{X — /i) = sin(A — /i), c(A — /i) = sm{h'y) 



(4.10) 



as that of the XXZ-spin chain. It is an easy check that ( f4.9|) and ( [4.8D thus found satisfy all 
YBE's ( ^^ and therefore, may be used for the braided QYBE's (^) and (|4l7| ) applicable 
to the periodic quantum mKdV model. Moreover due to the fact that i?-matrix ( ^4.9| ) commutes 
with Z- matrix ( [4.8| ), the braided QYBE's in the present case simplify further to 



and 



-Ri2(A - /i)Lij(A)L2j(Ai) = L2j{ii)Lij{X)Ri2{X - ^) 
i?i2(A - ii)T,{X)Z^2^TM = TMZ{2^T^{X)R^2{X - /i), 



(4.11) 
(4.12) 
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respectively. 

Before proceeding further we must justify our choice of the monodromy matrix taken in 
the same form as in the ultralocal case: 



Ta{X) = LaN{X) ■ ■ ■ ^a,fc+l(A)L„,fc(A) . . . L,i(A). (4.13) 

Evidently this is possible when the nonultralocal property does not spoil the detailed structure 
of this global object and one gets Lk^iLk{^) ~ LkLk+i{^). From ( |3.1|) we obtain the exphcit 
form 

Lk+iLk{i)=\ T{k+i,k) (4-14) 



^21 ^22 



where 



k 



4^''^ = -^^^m^^,r\w^)-' + wi:^,{w^r') (4.i5) 



Using nonultralocal relations ( |3.2| ) one shows that for all expressions ([4.15|) , L^-^^^''^'* = q L^j''''^^^ 



hold, which leads to the conclusion Lfc+iLfc(^) = q LkLk+i{^) justifying the choice ( |4.13| ). 



For identifying the commuting set of conserved quantities from ( |4.7| ) one may find the k- 
trace following the prescription of |21|, [T5|. However since in our case the i?- matrix commutes 
with Z we may start from ( [4. 121 ), multiply it by Z and take the trace as 



tr{Z^2Tl{\)Z{^T2{^l)) = tr{ZuT2{ij)Z^^'Ti{X)), 
forcing the standard relation 

[trT{\),trT{fi)] = . 

Therefore for 

r(A) = trT{X) = A{X) + D{X) becomes the generator of the conserved quantities with the 
Hamiltonian of the mkdV system produced as 5^'/nr(A) | (A = oo) = C3 = H at the continuum 
limit. Note that since the symmetry of the i?-matrix allows [Ri2,kik2] = with k = q~'^'^^ , 
one may also consider 

r(A) = tr{kT{X)) = q-^A{X) + q^D{X) 
also for generating the commuting family of conserved quantities. 
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5 Algebraic Bethe ansatz solution for the eigenvalue 
problem 

For solving the eigenvalue problem of the quantum mKdV model 

r(A) I Ai,A2,- ■■,A„^ > = {q-^ A{X) + q^D{X)) \ Xi,X2,- ■■ ,Xm> 

= A(A) I Ai,A2,---,A„. > (5.1) 

including that of its energy spectrum we may proceed in a way similar to the standard algebraic 
Bethe ansatz formalism , keeping at the same time the track of the nonultralocal signatures. 
Considering B{X) and C(A) as the creation and anihilation operators of the pseudoparticles 
we assume the eigenvectors in the usual Bethe ansatz form 

I Ai, A2, ■ ■ ■ , >= B{Xi)B{X2) . . . B{X^) I > (5.2) 

where | > is the pseudovacuum. For calculating the eigenvalue of r(A) = q~'^A{X) + q'^ D(X) 
we need to know first the commutation relations of ^(A) and D{X) with B[n), which may be 
obtained from the braided QYBE ( |4.12| ) using the R, Z-matrices given as ( |4.9| ) and ( |4.8[ ): 

D(X)B(^,) = ,2^|A^B(^)D(A) - ,=^|^BWB(A) (5.3) 

Note that the multiplicative factors q^"^ appearing in the above expressions is the contribu- 
tion from nonultralocality. For proceeding further along the QISM line one must define a 
pseudovacuum such that 

y..^ M(A) B{X)\ (a^{X) *\ 

- [ C(A) D(A) j ' ° ^" 1 P^{X) ) ' ° ^^-^^ 



which corresponds to the action of Lj on local | >j, where | >= ni=i I >i 

as Li \ >i= ( ^^^},.^ ] I >j . That is the pseudovacuum must be the eigenstate of 

the diagonal elements of Lj and anihilated by one of its off-diagonal elements. 

However as evident from the form of ( ^.l] ) naive application of a single Lax operator 
can not construct such a pseudovacuum. The situation is more like the Toda chain ||2^ or 
the lattice Liouville model [H[. Therefore one may either follow the functional Bethe ansatz 



method formulated in |^ or the remedy suggested in 0] . We adopt here the brilliant solution 
of by considering the action of a fused Lax operator L*^^^ = LfcL^+i in constructing the 
pseudovacuum. 

A careful analysis involving the action of the fused Lax operator elements ( [4.15 ) shows that 
here the pseudovacuum may be chosen in the form 

N 

1 >= n 1 >fc 

(^=1) 
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where 



I >,= - qW:)5{W,_,, + {W,)-'). (5.5) 

For checking this claim and deriving the required relation ( |5.4| ) we consider the action of 
^(A:+i,fe) ( |4.15| ) on ( |5.5| ). This gives for the diagonal elements 

= - (i - (Ao'g) - + (w^fc^)"') (5.6) 



I.e. 



a(A) = 22e^(^+^) sin(A + ^) (5.7) 



with = e*'^ and similarly 

Lg"-'^'^ 1 = + (Ae)2(w^,V)"'w^fc') 1 > 



yielding 



1-(A0V' |0>fc (5.^ 



= 22e*(^-^) sin(A - ^). (5.9) 
On the other hand for the off-diagonal element we obtain 

I > = -zAe [m+,)-\W,)-' + W^,+i(W^,+)-^) I 

= -zAe (g-^(M^,+)-^(W^.-)-^ - iW^)-\W^)-') I >,= (5.10) 

using the relation {W^)^^{W^)^^ = q^^ (Wj^)^^ {W,r)^^ , which corresponds to the anihilation 
operation. The other element however gives l[2^^''^'* | >k^ 0. The eigenvalue problem (|5.1| ) 



of the mKdV model therefore may be solved using the relations ( p.3|) taking into account the 
pseudovacuum properties ( |5.4| ). Finally inserting the results (|5.7| ) and (|5.9|) one obtains 

A(A) ^ ,-Hxfni^r'^^).,mxfii(r^) 

2 "-Ji sm(A - A^) 

+ ,-(-l^f-^-)sm"(A-;^)n "°<.\-/'tf' 1 (5-11) 

2 j^-^ sm(A — Ajj 

The determining equation for the rapidity parameters Xj may be obtained from the vanishing 
residue condition of the eigenvalue equation ( ^.11| ) in the usual way: 

-.-,N-.m / sin(A+f) \ _ - / sin(A-A, + 7n) \ ..... 
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Note that the nonultralocal property of the model is reflected in the q~^^ factors appearing in 
the above relations. Such factors were found to appear also in ultralocal models with twisted 
boundary conditions |23, 24]. However for the mKdV model considered here, they depend 
on the total number of lattice points 2N as in 0] and more interestingly, on the number of 
pseudoparticle excitations m. Comparison of ( ^.12|) with the Bethe equation for the spin-| 
XXZ chain with twisted boundary condition 

/ sm(A + ^) \ " _ " / sm(A-Aj + 7S) \ 

reveals the startling fact that, if the extra q factors are ignored the mKdV model on a lattice 
with 2N points becomes equivalent to the spin-| XXZ chain with total lattice points A^. Of 
course the hamiltonian of the respective systems corresponds to different conserved quantities 
and their expansion point of spectral parameters for generating conservation laws also differs. 
For obtaining the mKdV field model one has to consider ultimately the continuum limit. 



6 Concluding remarks 

The spin- 1 XXZ chain is known to be connected with the six- vertex model, a classical sta- 
tistical model in 2-dimrnsions. Therefore since the finite-size corrections and the conformal 
properties of the six-vertex model is a well studied subject, the intriguing relation between 
the quantum mKdV and the spin-| XXZ chain found here should be an important result in 
analyzing the conformal properties of the mKdV model with q = . Let us recall that the 
six-vertex model (with a seam) yields p4| for the ground state energy 



and for the excited states 



E^^ -E, = |(A + A) + O(^) - Po = f (A - A) + 0{^) 

where c = 1 t^ttt, ^ = 2,3, .. . for k = 1 is the central charge and A, A are conformal 

weights of unitary minimal models of the corresponding conformal field theory. 

The conformal properties of the mKdV model have enhanced importance due to its mapping 
into the KdV model through Miura transformation. The KdV field in turn is related to the 
generators of the Virasoro algebra answering to the 2-d conformal transformations [^, ^ 



In this context the continuum limit of the fused Lax operator ( [4.14| ), which plays a crucial 



role in constructing the pseudovacuum and deriving the Bethe equation ( |5.12| ), is also worth 
analyzing. The linear system corresponding to this Lax operator is $^+2(0 — L^'^\0^kiO > 
which at the continuum limit would yield the relations 



^ / -/ N 3n ^ /„N , ^ /„N / U[X) 



MO) = -^iv-{x)a')m, and $..(0) = u(x) ) ^^^^ ^^'^^ 
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where $(0) = = 0) and the fields u, u are given by 

u{x) = iv~ [x] + v~ (x), u{x) = —iv~ {x) + v~ {x) (6.2) 

i.e. they represent KdV fields related trough Miura transformation. Replacement of v~ ^ 
would generate similarly KdV fields through . Due to such relations between mKdV and 
KdV fields and their connection with the Liouville field the recently discovered equivalence 
between quantum lattice Liouville model and spin (— |) XXZ chain acquires renewed 
importance. 

It is also hoped that the main results of the present work, e.g. construction of the exact 
lattice Lax operator for the quantum mKdV model, the R and Z matrices as solutions of the 
braided QYBE, the Bethe ansatz solution of the related eigenvalue problem and finally the 
mapping between quantum mKdV and spin-| XXZ chain, would be of significance in the 
alternate approach of ||26| in formulating CFT through quantum KdV related theory. 
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